Two-photon path-entangled states in multi-mode waveguides 
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We experimentally show that two-photon path-entangled states can be coherently manipulated 
by multi-mode interference in multi-mode waveguides. By measuring the output two-photon spatial 
correlation function versus the phase of the input state, we show that multi-mode waveguides perform 
as nearly-ideal multi-port beam splitters at the quantum level, creating a large variety of entangled 
and separable multi-path two-photon states. 
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Quantum states of photons distributed between sev- 
eral optical paths play the central role in linear opti- 
cal quantum computation In particular, entangled 
states, where no decomposition into a product of cither 
single-path or single-photon states exists, are of great 
importance due to their non-local nature 0, Q • Entan- 
gled photons states are generated either by non-linear 
processes in crystals by cascaded emission of single 
quantum emitters Q , or by quantum interference of non- 
classical light on a beam splitter @, 0] ■ 

Quantum interference in arrays of beam-splitters (BSs) 
is used in many linear optical quantum computation 
schemes to create and manipulate multi-photon multi- 
path states [8 10]. The simplest example is the Hong- 



Ill. In this fundamcn- 



Ou-Mandcl (HOM) interference 
tal effect, two identical photons input each on a different 
port of a 50:50 BS, cannot exit in two different ports 
due to destructive interference of the two possible paths. 
Therefore, they always exit bunched together on either 
one of the output ports, in an entangled state. A ma- 
jor obstacle on the way to the implementation of quan- 
tum optical circuits is the large number of BSs required, 
and the increasing complication of their alignment. One 
way to overcome this problem is through miniaturiza- 
tion of the optical circuitcs. Indeed, quantum interfer- 
ence has been recently demonstrated in integrated op- 
tical circuites composed of evanescently coupled single- 
mode waveguides embedded in solids A con- 
ceptually different route towards robust implementation 
of quantum optical circuites may come in the form of 
multi-mode interference (MMI) devices [Hj]. These com- 
pact replacements for BS arrays, usually based on planar 
multi-mode waveguides (MMWs), arc already used ex- 
tensively in modern classical optical communication net- 
works. They have also been proposed to be useful for cre- 
ation and detection of multi-photon states [l(|, as they 
naturally implement Bell multiport BSs Q. A step to- 
wards their use in quantum networks was recently made 
with the demonstration of HOM interference of two sep- 
arated photons in an MMW [la ]. The question still re- 
mains, however, whether entangled states, which unlike 
separated photons, carry relative phase information, can 
be coherently manipulated by MMWs. 



Here we experimentally show that the answer to this 
question is positive. We utilize MMI in a two-mirror, tun- 
able planar MMW [l?! to implement multi-port BS ar- 
rays of up to 5 input and 5 output ports, and explore the 
propagation of non-classical, path-entangled two-photon 
states through them. We measure the two-photon cor- 
relations at the output, and find that, for any relative 
phase between the two paths, they agree very well with 
the multi-path, two-photon states expected at the output 
of ideal multi-port BSs. 

The experimental system is described in Fig. HJa). 
States of the form ^ = -±= (|2, 0) + e l *|0, 2)) are created 
by first splitting a 404 nm continuous-wave (cw) diode- 
laser beam into two arms. The relative phase between 
the beams, (f>, is controlled and stabilized by a piezoelec- 
trically movable mirror on one arm [see Fig. QJa)]. The 
state in the two arms at this point is \a, e^a), where 
a represents a coherent state of average photon num- 
ber |a| 2 . Note that the single photon part of this state 
is tpf = (| 1, 0> + e l *|0, 1)). The two beams then un- 
dergo type I collinear degenerate spontaneous parametric 
down-conversion (SPDC) in a properly oriented 2 mm 
long /3-Barium-Boratc (BBO) crystal. In this process, 
each 404 nm photon transforms into two 808 nm photons 
propagating together. Since SPDC is much faster than 
the optical period, the relative phase between the two 
paths is kept as it was before the BBO crystal. Spectral 
filters are used to eliminate all remaining photons out- 
side a 3 nm band centered at 808 nm. The two photon 
part of the state after the BBO and the spectral filters 
is thus ip2 ■ There is no single-photon part, and for low 
enough pump powers, the content of the next order - the 
four photon part, is negligible. 

The beams are then inserted into an MMW made of 
two parallel metallic mirrors 17|. The distribution of 
one component of the optical field between the mirrors, 
assuming perfect reflection and paraxial propagation, is 
given by [15j , 



E(x, z) = e- lkz An sin [mr (x - D/2) / D] e i2 ™V^ 
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FIG. 1: (a) The experimental setup. DL, 404 nm cw Diode 
laser; SF, Spatial filter; HWP, Half wave plate; PBS, Polar- 
izing beam splitter; RR, Retroreflector; QWP, Quarter wave 
plate; M, Mirror; BBO, /?-Barium-Borate crystal; BPF, Band- 
pass spectral filter; MMW, Multimode waveguide; BS, beam 
splitter; MMF, Multimode fiber; APD, Avalanche photodi- 
ode; CE, Correlation electronics. The thin purple (wide red) 
lines represent the pump (down conversion) beams. Movable 
parts are marked by black double arrows, (b) Calculated in- 
tensity distribution in an MMW for an input beam centered 
about x = D/A. The white boxes mark propagation lengths 
where various 2x2 BSs are implemented, (c)-(f) Schematic di- 
agrams showing the configurations of input and output ports 
used in the experiments. 



with k = 2ir/\ and zq = 8D 2 /X, where D is the width of 
the MMW, A is the wavelength of the incident light, and 
the constants A n are determined by the incident field dis- 
tribution. From Eq.[T]it can be shown that for any input, 
full imaging occurs at z = zq, reflection about the cen- 
ter of the MMW occurs at z = Zo/2, and two-way, equal 
beam splitting with a relative phase of7r/2 (— it/ 2) occurs 
at z = zo/4 (z = 3zp/4). These effects are known as gen- 
eral interference [15(. Unequal BSs, and BSs of more than 
two ports are also possible, but only for specific positions 
of a localized input beam. This is known as restricted in- 
terference [lij]. For example, for an input beam localized 
about x = zLD/4, an unequal two-way BS will be real- 
ized at z = (2m — l)zo/8, where m is an integer. This 
is readily seen in Fig. (Hb), where the intensity distribu- 
tion inside the MMW, calculated for a Gaussian input 
beam localized about x = D/4, is presented. Generally, 



an NxN BS would be generated by restricted interfer- 
ence for beams localized about x — (2p — 1 — N)D/2N, 
for z — qzq/AN, where p and q are integers, and p goes 
from 1 to N. For an even q, the BS is equal. Using 

Eq. [Tj the transition matrix, = ^T^J , can be 

calculated for any N This matrix applies for co- 

herent beams as well as for a single photon. In order 
to obtain the transition matrices for states of M identi- 
cal photons, one should sum up all the products of the 
single-photon matrix elements relevant to the transition 
between an initial, V, and a final, fl, Af-photon configu- 
ration, 

i~\t~ m 

(<0^=v§£n <) ■ < 2 > 

where represent permutations of the initial configu- 
ration v, and N$ (Nn) is the number of all different 
permutations of the initial (final) configuration. The 
calculation of these iV-port, Af-photon transition ma- 
trix elements is, in general, a computationally difficult 
task, equivalent to the calculation of the permanent of 
a matrix. However, it is feasible for small enough Af 
and N, and for N = 2 there even exists an analytic ex- 
pression for any Af [l9|. Given an initial wavefunction 
(in terms of the initial amplitudes for each distribution 
of the Af photons among the N input ports), the final 
state, ipout, is readily found, and from it, any correlation 
function can be calculated. In particular, for a general 
two-photon state, the two-photon correlation probability 
between ports m and n is just 

P$ n = \( m ,n\ip out )\ 2 , (3) 

where \m, n) is a state where one photon is on the m th 
output port, and the other is on the n th output port. 

In our experimental setup, the length of the MMW 
is fixed by the length of the mirrors, L. However, the 
adjustment of the relative length, £ = L/zq = LX/8D 2 , 
is possible through the adjustment of the separation be- 
tween the mirrors. For the experiments presented below 
the MMW was adjusted to function as 2x2, 3x3, 4x4, 
and 5x5 BSs. Figs. [IJc)-(f ) show the corresponding con- 
figurations of input and output ports used in these ex- 
periments. We note that due to the finite conduction 
of the mirrors the MMW has a slightly different relative 
length for different polarizations. To avoid any possible 
complications that may arise, the polarization of the in- 
cident photons (in both beams) was set parallel to the 
mirrors, such that only TE modes were excited. In or- 
der to measure two-photon correlations at the output of 
the MMW, the light at its output facet is split and im- 
aged on two multi-mode fibers, each connected to a single 
photon detector (Si avalanche photodiode). The digital 
signals from the two detectors are input into correlation 
electronics that shorten and multiply them, yielding the 
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rate of coincidence events in a 7 ns time window [20j • As 
the fibers can be placed in front of any pair of output 
beams [see Fig. QJa)], both auto- and cross-correlations 
can be measured. Since in this method, the beams are 
split into two also when cross-correlations are measured, 
the measured rate of cross-correlations would be half that 
predicted using Eq. [3l For the comparison of theory and 
experiment we therefore define the modified correlation 
probability, 



2) 
m.n 



p(2) 
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.»/(2 -Sm,n). (4) 

Fig. [2la) presents the calculated modified two-photon 
correlation probabilities (Eq. 2]) versus the phase <f>, for 
three different relative propagation lengths at which 2x2 
BSs are implemented [see Fig. [Ub) and (c)]: C, = 1/4 
- a 50:50 BS with a relative phase of 7r/2, £ = 3/8 - a 
cos 2 7r/8 : sin 2 7r/8 (~85:15) BS with a relative phase of 
7r/2, and £ = 1/2 - a reflection about the center of the 
MMW (a '0:100 BS')- Note that all three cases are a part 
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of the same series of transformations given by (jl 

where \ is a 15:85 BS implemented for £ = 1 /8, and 
that for q = 8 (£ = 1), full imaging is obtained [13, [H|. 
The state ip® ^ s invariant only under reflection, and thus 
returns to its initial form only when full reflection occurs 
- once every four applications of T^l ( e -§- a ^ C=l/2)- 
Note that at the middle of this period (on the 50:50 BS at 
£=1/4) this state transforms into two separated photons, 
exhibiting inverse HOM interference. In contrast, the 
state ip2 i s unchanged (up to a global phase) for each 
and every application of T%\ ■ This is because this state is 
invariant under any unitary two-way beam splitter with 



a relative phase of ±7r/2 [21 



Fig. HJb) shows the measured correlation functions, 

(21 

g-mn (symbols), defined as the rate of coincidence events 
between output ports m and n, normalized by the ex- 
pected accidental coincidence rate, r a = r m r n w, r m be- 
ing the single-count rate in port m, and w = 7 ns the 
coincidence measurement time window. The total inte- 
gration time is 7 sec for each data point. The lines are the 
best fitted 2-7r-period sinusoidal functions. The measured 
background after the MMW (right panel of Fig. EJ[b), 
empty circles) is ~10 times larger than the expected ac- 
cidental coincidence rate. By blocking either one of the 
input beams, we have verified that indeed ^90% of this 
background is due to a small overlap between the col- 
lecting fiber and the neighboring output port. This is 
also confirmed by the increase of the background level 
when the beams are brought closer (see, e.g., Fig. HJb) 
below). The bare visibility of the oscillations measured 
for the 50:50 BS configuration (Fig. HJb), left panel) is 
72±3%, violating the classical bound of 50% [22j by more 
than 7 standard deviations. This indicates that the light 
after the MMW is still entangled, even without back- 
ground reduction. With the reduction of the measured 
background, the visibility reaches 83±3%. 




• 1-1 

■ 2-2 

o 1-2 

O OOOO OOP oo o<j > 
0.5 1 0.5 1 0.5 1 

0/2JI 



=0 100 



=Jt 100 




FIG. 2: (a) Two-photon correlations vs. the phase, <f>, calcu- 
lated using Eq. Ufor relative waveguide lengths of 1/4 (left), 
3/8 (center), and 1/2 (right), where the MMW functions as 
different types of 2x2 BSs, as illustrated in Fig. [He). The 
correlation between ports m and n is marked by 'm-n'. (b) 
Measured correlation functions. The symbols are the mea- 
sured values and the lines are best fitted 2n period sinusoidal 
functions, (c) Corresponding two-photon correlation maps 
extracted from the fitted curves for (f> — (top) and (f> — n 
(bottom). 



Fig. [He) presents 'correlation maps' for each of the 

(2) 

relative lengths, showing the elements g m h for two initial 
phases, <f> = (top) and <f> = 7r (bottom), as extracted 
from the fits to the measured data. For each phase, the 
expected periodicity with respect to the relative length 
of the MMW is clearly seen. 

We proceed to examine 3x3 BSs. The input and out- 
put ports used are as illustrated in Fig.QJd). Fig.[3]shows 
the two possible cases. On the left we present the calcu- 
lated (a) and measured (b) correlations of a 3x3 equal 
BS (£ = 1/3), while on the right the correlations for an 
unequal BS (£ = 5/12) are presented. For the equal BS, 
the six possible two-port correlation measurements di- 
vide into three pairs, oscillating with a phase difference 
of 27r/3 between them. Each pair contains an autocorre- 
lation of one port and the cross correlation of the other 
two ports. There are six special phases where two such 
correlations pairs meet. For three of them there is an 
enhancement of one pair over the other two, while for 
the other three, one pair is completely suppressed. This 
is further visualized in Fig. 03c), where the correlation 
maps for </> = and (f> = it are shown. In contrast to 
the equal 3x3 BS, the phase dependencies of the corre- 
lations in the unequal 3x3 BS show only two different 
phases of oscillation and only partial visibilities. How- 
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FIG. 3: (a) Calculated two-photon correlations vs. the phase, 
cj>, for two cases of 3x3 multi-port BSs implemented by set- 
ting the relative length of the MMW to 1/3 (left), and 5/12 
(right). The locations of the input and output ports are il- 
lustrated in Fig. [Hd). The port configurations ('m-n') of all 
correlations described by a certain curve are listed next to it. 
(b) The measured correlations (symbols) and their sinusoidal 
fits (lines), (c) Correlation maps for — and = m. 



ever, the two cases can be related one to the other if one 
notes that the unequal BS is composed of an equal 2x2 
BS (( = 1/4) on the two outer ports, followed by an equal 
3x3 BS of C = 1/6, different from that of C = 1/3 only in 
the order of the relative phases [l8|. This is best seen in 
the correlation maps shown in Fig. |3{c). For = 0, the 
cross-correlation between the two outer arms switch val- 
ues with the corresponding auto-correlations, while for 
= 7r there is no difference between the two 3x3 BSs, 
due to the invariance of -02 on 2x2 ir/2 BSs. 

To examine the effect of even more complex BSs on the 
states tjj^, we performed measurements for equal 4x4 and 
5x5 BSs. The input and output ports used are as illus- 
trated in Figs. [IJc) and (f), respectively. Fig. HJa) [(b)] 
presents the calculated [measured] phase dependence of 
these correlations, and Fig. @|c) shows the correspond- 
ing correlation maps for = and = 7r. In the case 
of the 4x4 equal BS, all the auto correlations and the 
cross-correlations of symmetric ports oscillate together 
with a phase of 0, and all the cross-correlations of asym- 
metric ports oscillate together with a phase of tt. The 
state created for = therefore contains no asymmetric 
cross correlations, while that created for = ir contains 
only asymmetric cross correlations. This is visualized in 
the correlation maps of Fig. B|c). For the equal 5x5 
BS, analogously to the equal 3x3 BS, the 15 correlations 
group in 5 triplets that oscillate with phase differences 



FIG. 4: (a) Calculated two-photon correlations vs. the phase, 
0, for an equal 4x4 multi-port BS (left), implemented by set- 
ting the relative length of the MMW to 1/8, and an equal 5x5 
multi-port BS (right), for which the relative length is 1/5. 
The locations of the input and output ports are illustrated in 
Fig. [He) and (f), respectively. The output ports configura- 
tions ('m-n') of all correlations described by a certain curve 
are listed next to it. (b) Measured correlation functions. Full 
black symbols represent auto correlations, while gray-filled 
and empty symbols represent cross-correlations. The solid 
lines are sinusoidal fits, (c) Corresponding two-photon corre- 
lation maps for <j> = and = tt. 



of 27r/5. Here too, some correlations can be either rela- 
tively enhanced or completely suppressed, as visualized 
in Fig. gjc). 

In summary, we experimentally show, for the first time, 
that path-entangled quantum light can be coherently ma- 
nipulated by multi-port BSs, naturally implemented by 
MMI in MMWs. The manipulation is shown to maintain 
the information on the relative phase between the two 
paths. This persists even for complex manipulations with 
a high number of output ports, allowing for the robust 
creation of a large space of two-photon multi-path states, 
controlled by the initial phase. MMI is thus shown to be 
a robust and simple approach for the implementation of 
various quantum optical circuites of high complexity. 
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